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CN 1 ABSTRACT 

Using a directional spherical wavelet analysis we detect the integrated Sachs-Wolfe (ISW) 
effect, indicated by a positive correlation between the first-year Wilkinson Microwave 
<N ; Anisotropy Probe (WMAP) and NRAO VLA Sky Survey (NVSS) data. Detections are made 

using both a directional extension of the spherical Mexican hat wavelet and the spherical but- 
terfly wavelet. We examine the possibility of foreground contamination and systematics in the 
WMAP data and conclude that these factors are not responsible for the signal that we detect. 
£vq . The wavelet analysis inherently enables us to localise on the sky those regions that contribute 

■ most strongly to the correlation. On removing these localised regions the correlation that we 

\Q \ detect is reduced in significance, as expected, but it is not eliminated, suggesting that these 

regions are not the sole source of correlation between the data. This finding is consistent 
with predictions made using the ISW effect, where one would expect weak correlations over 
the entire sky. In a flat universe the detection of the ISW effect provides direct and indepen- 
dent evidence for dark energy. We use our detection to constrain dark energy parameters by 
O \ deriving a theoretical prediction for the directional wavelet covariance statistic for a given 

cosmological model. Comparing these predictions with the data we place constraints on the 
^ ■ equation-of- state parameter w and the vacuum energy density Qa. We also consider the case 

of a pure cosmological constant, i.e. w = -1. For this case we rule out a zero cosmological 
constant at greater than the 99.9% significance level. All parameter estimates that we obtain 
are consistent with the standand cosmological concordance model values. Although wavelets 
rN \ perform very well when attempting to detect the ISW effect since one may probe only the 

^ ■ regions where the signal is present, once all information is incorporated when computing pa- 

rameter estimates, the performance of the wavelet analysis is comparable to other methods, as 
expected for a linear approach. 

Key words: cosmic microwave background - cosmology: observations - methods: data anal- 
ysis - methods: numerical. 
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1 INTRODUCTION 

Strong observational evidence now exists in support of the A cold 
dark matter (ACDM) fiducial model of the universe; specifically, a 
universe that is (nearly) flat and dominated by an exotic dark en- 
ergy component. We know very little about the origin and nature 
of this dark energy, but we now have strong evidence in support of 
its existence and relative abundance. Much of this evidence comes 
from recent measurements of the cosmic microwave background 
(CMB) anisotropics of high resolution and precision, in particular 
the recent Wilkinson Microwave Anisotropy Probe (WMAP) data 
(Bennett et al. 2003a). The existence of dark energy has also been 
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independently found by measurements of the luminosity distance 
to Supernova Type la (Riess et al. 1998; Perlmutter et al. 1999). 

At this point, the confirmation of the fiducial ACDM model 
and the existence of dark energy by independent physical methods 
is of particular interest. One such approach is through the detection 
of the integrated Sachs-Wolfe (ISW) effect (Sachs & Wolfe 1967). 
CMB photons are blue and red shifted as they fall into and out 
of gravitational potential wells, respectively, as they travel towards 
us from the surface of last scattering. If the gravitational poten- 
tial evolves during the photon propagation, then the blue and red 
shifts do not cancel exactly and a net change in the photon energy 
occurs. This secondary induced CMB anisotropy (the ISW effect) 
exists only in the presence of spatial curvature or, in a flat universe, 
in the presence of dark energy (Peebles & Ratra 2003). The recent 
WMAP 1-year data has imposed strong constraints on the flatness 
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of the universe (Bennett et al. 2003a; Spergel et al. 2003), hence 
any ISW signal may be interpreted directly as a signature of dark 
energy. 

It is difficult to separate directly the contribution of the ISW 
effect from the CMB anisotropics, hence it is not feasible to de- 
tect the ISW effect solely from the CMB. Instead, as first proposed 
by Crittenden & Turok (1996), the ISW effect may be detected by 
cross-correlating the CMB anisotropics with tracers of the local 
matter distribution, such as the nearby galaxy density distribution. 
A detection of large-scale positive correlations is a direct indication 
of the ISW effect, and correspondingly, evidence for dark energy. 

The first attempt to detect the cross-correlation be- 
tween the CMB and the nearby galaxy distribution was 
performed by Boughn & Crittenden (2002) using the Cos- 
mic Background Explorer-Differential Microwave Radiometer 
(COBE-DMR; Bennett et al. 1996) and NRAO VLA Sky Survey 
(NVSS; Condon et al. 1998) data. 1 No cross-correlation was found; 
Boughn & Crittenden (2002) conclude that a future experiment 
with better sensitivity and resolution was required to make any 
detection. The WMAP mission has now provided a suitable ex- 
periment, and several groups have since reported detections of the 
ISW effect at a range of significance levels using various tracers 
of the local universe. Boughn & Crittenden (2004, 2005) cross- 
correlate the WMAP data with two different tracers of the nearby 
universe: the hard X-ray data provided by the High Energy Astron- 
omy Observatory- 1 satellite (HEAO-1; Boldt 1987) and the NVSS 
data. They make a statistically significant detection of the cross- 
correlation at the 1.8-2.8cr level at scales below 3°. Nolta et al. 
(2004) perform an independent analysis of the WMAP and NVSS 
data and confirm the existence of dark energy at the 95% sig- 
nificance level. Fosalba & Gaztanaga (2004) cross-correlate the 
WMAP data with the APM galaxy survey (Maddox et al. 1990) 
and report a cross-correlation detection at the 2.5cr level on scales 
of 4°-10°. Fosalba etal. (2003) cross-correlate the WMAP data 
with the Sloan Digital Sky Survey (SDSS; York etal. 2000) and 
detect an ISW signal at the 3<r level. These same two data sets 
are also analysed by Scranton et al. (2003). 2 Afshordi et al. (2004) 
cross-correlate the WMAP data with the near infrared Two Mi- 
cron All Sky Survey Extended Source Catalogue (2MASS XSC; 
Jarrett et al. 2000) and detect an ISW signal at the 2.5<x level. 
Padmanabhan et al. (2004) make a 2.5<x detection of the ISW ef- 
fect by cross-correlating the WMAP and SDSS data, and subse- 
quently use the detection to constrain cosmological parameters. 
Other works have focused on the theoretical detectability of the 
ISW effect for various experiments and, in some cases, the use 
of such detections to constrain cosmological parameters (Afshordi 
2004; Hu & Scranton 2004; Pogosian 2004; Pogosian et al. 2005; 
Corasaniti et al. 2005). 

The previous works discussed all perform the cross- 
correlation of the CMB with various tracers of the near uni- 
verse large scale structure (LSS) in either real or harmonic space, 
using the real space angular correlation function or the cross- 
angular power spectrum respectively. Recently, Vielva et al. (2006) 
adopt a different measure by performing the cross-correlation in 
spherical wavelet space using the azimuthally symmetric spher- 



As the detection of the large-scale effect is cosmic variance limited, one 
requires (near) full-sky maps. 

2 See Afshordi et al. (2004) for a critical discussion of the analyses done 
by Fosalba & Gaztanaga (2004), Fosalba et al. (2003) and Scranton et al. 
(2003). 



ical Mexican hat wavelet (SMHW). Spherical wavelets have al- 
ready been used in many astrophysical and cosmological applica- 
tions. For example, spherical wavelets have been used extensively 
to test the CMB for non-Gaussianity and isotropy (Barreiro et al. 
2000; Cayonetal. 2001, 2003, 2005; Martmez-Gonzalez et al. 
2002; Vielva etal. 2004; Mukherjee & Wang 2004; Cruz et al. 
2005, 2006; McEwen et al. 2005b; Liu & Zhang 2005; Wiaux et al. 
2006). However, Vielva et al. (2006) present the first use of spher- 
ical wavelets for cross-correlating the CMB with tracers of LSS 
in an attempt to detect the ISW effect. Since the ISW effect is lo- 
calised to certain (large) scales on the sky, wavelets are an ideal 
tool for searching for cross-correlations due to the inherent scale 
and spatial localisation afforded by a wavelet analysis. Vielva et al. 
(2006) examine the covariance of the SMHW coefficients of the 
WMAP and NVSS data for a range of scales, making a detection 
of the ISW effect at the 3.3<x level at scales on the sky of 6°-8°. 
Moreover, the detection is used to constrain cosmological parame- 
ters that describe the dark energy. 

There is no physical reason to assume that the local corre- 
lated structures induced in the CMB anisotropics by the near LSS 
are rotationally invariant; indeed, it is known that Gaussian random 
fields are characterised by features that are not necessarily rota- 
tionally invariant Barreiro et al. (1997, 2001). Thus, other spheri- 
cal wavelets that are not azimuthally symmetric, i.e. directional 3 
wavelets, may be equally (or more) suitable for probing the data 
for cross-correlations. Herein we extend the analysis performed 
by Vielva et al. (2006), using directional wavelets to examine the 
cross-correlation of the WMAP and NVSS data. We use any detec- 
tions of the ISW made to constrain cosmological parameters that 
describe the dark energy. 

The remainder of this paper is structured as follows. The direc- 
tional continuous spherical wavelet transform and the wavelets that 
we consider are described in section 2. In section 3 we define the 
wavelet covariance estimator used to test for correlations and also 
define the theoretical covariance predicted for a given cosmological 
model (the derivation for which is presented in Appendix A). We 
then compare the expected performance of various wavelets for de- 
tecting the ISW effect. In section 4 we give a brief overview of the 
data considered and the analysis procedure. Results are presented 
and discussed in section 5. We describe the detections made and 
then use the detections to place constraints on dark energy parame- 
ters. Concluding remarks are made in section 6 



2 CONTINUOUS SPHERICAL WAVELET TRANSFORM 

To perform a wavelet analysis of full- sky maps defined on the 
celestial sphere, Euclidean wavelet analysis must be extended to 
spherical geometry. A wavelet transform on the sphere has al- 
ready been constructed and developed (Antoine & Vandergheynst 
1998, 1999; Antoine et al. 2002, 2004). 4 The spherical wavelet 
construction is derived entirely from group theoretic principles, 
however recently Wiaux et al. (2005) reintroduce the formalism 
in an equivalent, practical and self-consistent approach that is in- 
dependent of the original group theoretic framework. We adopt 
this latter approach and extend the decomposition to anisotropic 

3 We use the term 'directional wavelet' to refer to a wavelet that is not 
azimuthally symmetric. This may differ to the definition of a directional 
wavelet used by some other authors. 

4 This framework has been extended to wavelet frames on the sphere also 

by Bogdanova et al. (2005). 



Detection of the ISW effect and corresponding dark energy constraints 3 



dilations (McEwen et al. 2005a). Moreover, we also apply a fast 
algorithm for performing the wavelet transform on the sphere 
(McEwen et al. 2005a), based on the fast spherical convolution al- 
gorithm developed by Wandelt & Gorski (2001). We present here 
a brief overview of the continuous spherical wavelet transform 
(CSWT), but refer the reader to our recent work (McEwen et al. 
2005a) for more details on the analysis and fast algorithms. 



2.1 Wavelet transform 

The correspondence principle between spherical and Euclidean 
wavelets developed by Wiaux et al. (2005), relates the concepts of 
planar Euclidean wavelets to spherical wavelets through a stereo- 
graphic projection. The stereographic projection is used to define 
affine transformations on the sphere that facilitate the construction 
of a wavelet basis on the sphere. The spherical wavelet transform 
may then be defined as the projection on to this basis, where the 
spherical wavelets must satisfy the appropriate admissibility crite- 
rion to ensure perfect reconstruction. 

The stereographic projection is defined by projecting a point 
on the unit sphere to a point on the tangent plane at the north pole, 
by casting a ray though the point and the south pole. The point on 
the sphere is mapped on to the intersection of this ray and the tan- 
gent plane (see Fig. 1). The stereographic projection is radial and 
conformal (Wiaux et al. 2005), hence local angles are preserved un- 
der the transform, i.e. directionality is preserved. The stereographic 
projection operator that preserves the L 2 -norm of functions is de- 
noted II, with inverse II" 1 . 

The stereographic projection operator may be used to ex- 
tend the concept of dilations to the sphere. In particular, we adopt 
here the extension of the isotropic dilation defined by Wiaux et al. 
(2005) to anisotropic dilations (McEwen et al. 2005a). 5 Dilations 
on the sphere are constructed by first projecting the sphere on to 
the plane using the stereographic projection, performing the usual 
Euclidean dilation in the plane, before re-projecting back on to the 
sphere using the inverse stereographic projections. The spherical 
dilation is thus defined by 



D(a,b) ^IT 1 d(a,b)U., 



(1) 



where d(a, b) is the anisotropic Euclidean dilation on the plane, de- 
fined for a function on the plane p e L 2 (R 2 ) by [d(a, b)p] (x,y) = 
a~ ll2 b~ ll2 p{a~ l x,b~ l y), for the non-zero positive scales a, b e R+. 
The L 2 -norm of functions is preserved by the spherical dilation op- 
erator since both the stereographic projection and the Euclidean 
dilation preserve the norm. We therefore obtain the following def- 
inition of the spherical dilation of a square-integrable function on 
the sphere s e L 2 (S 2 ): 



(2) 



[D(a,b)s](aj) = [A(a,b,6,(f>)Y /z s(aj w/b ) , 
where a) aJb = (6aj,,(/>ajb), 

tsm(0ajb/2) = tan(0/2) ^a 2 cos 2 + b 2 sin 2 

and tan(0 fl ^) = - a tan(0). The spherical coordinates with colatitude 
6 and longitude are denoted by co = (#, 0) e S 2 . The A{a, b, #, 0) 
cocycle term follows from the various factors introduced to pre- 
serve the L 2 -norm of functions. The cocycle of an anisotropic 
spherical dilation is defined by 




South pole 

Figure 1. Stereographic projection of the sphere onto the plane. 
4a 3 b 3 

A(a, b, 6, 0) = , (3) 

(A_ cos 6 + A + ) 2 

where 

A+ = a 2 b 2 ± a 2 sin 2 + b 2 cos 2 . 

Finally, we require an extension of Euclidean translations to 
the sphere in order to construct a wavelet basis on the sphere. 
The natural extension of translations to the sphere are rotations. 
These are characterised by the elements of the rotation group 
SO(3), which we parameterise in terms of the three Euler angles 
p = (a,j3, y). 6 The rotation of a function on the sphere is defined by 



[R(p)s](oj) = s{p- l aj\ p e SO(3) . 



(4) 



A wavelet basis on the sphere may now be constructed 
from rotations and dilations of a mother spherical wavelet 
G L 2 (S 2 ). The corresponding wavelet family on the sphere 
{^a,b,p ~ R(p)D(a,b)ifr, p e SO(3), a,b e R+} provides an over- 
complete set of functions in L 2 (S 2 ). The CSWT of a function on 
the sphere is given by the projection on to each wavelet basis func- 
tion in the usual manner, 



Js 2 



(5) 



where the * denotes complex conjugation and d£l(aj) = sin#d#d0 
is the usual rotation invariant measure on the sphere. In this work 
we do not consider the synthesis of a function on the sphere from its 
wavelet coefficients. Indeed, it is only possible to synthesise a func- 
tion from its spherical wavelet coefficients for the case of isotropic 
dilations (McEwen et al. 2005a). 

The transform is general in the sense that all orientations in 
the rotation group SO(3) are considered, thus directional structure 
is naturally incorporated. It is important to note, however, that only 
local directions make any sense on S 2 . There is no global way of 
defining directions on the sphere - there will always be some sin- 
gular point where the definition fails. 

2.2 Mother spherical wavelets 

Admissible mother spherical wavelets may be constructed from the 
stereographic projection of admissible mother Euclidean wavelets 
on the plane: 



5 A similar anisotropic dilation operator on the sphere has also been inde- 
pendently proposed by Tosic et al. (2005). 



6 We adopt the zyz Euler convention corresponding to the rotation of a 
physical body in a fixed co-ordinate system about the z, y and z axes by 
y, and a respectively. 
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^a)) = \nr 1 M((o). (6) 

Directional spherical wavelets may be naturally constructed in this 
setting - they are simply the projection of directional Euclidean 
planar wavelets on to the sphere. 

We consider three spherical wavelets in our subsequent cross- 
correlation analysis: the spherical Mexican hat wavelet (SMHW); 
the spherical butterfly wavelet (SBW); and the spherical real Morlet 
wavelet (SMW). These spherical wavelets are illustrated in Fig. 2. 
Each spherical wavelet is constructed by the stereographic projec- 
tion of the corresponding Euclidean wavelet onto the sphere, where 
the Euclidean planar wavelets are defined by 

and 

^r(x;k) = Re(e*' x/V5 e- ||x||2/2 ) 

respectively, where k is the wave vector of the SMW. The SMHW 
is proportional to the Laplacian of a Gaussian, whereas the SBW 
is proportional to the first partial derivative of a Gaussian in the x- 
direction. The SMW is a Gaussian modulated sinusoid, or Gabor 
wavelet. 

The dilation parameter of each wavelet may be related to an 
effective size on the size for the wavelet. For directional wavelets 
we define two effective sizes on the sky defined in orthogonal di- 
rections. The SMHW and SBW are both derived from a parent 
Gaussian function and thus for each direction have the same ef- 
fective size on the sky, defined by ^ SMHW SBW ( C .) ~ 2 V2 c u where 
i = {1, 2} and c t = {a, H=i,2 (McEwen et al. 2005b). The first effec- 
tive size of the SMW defines the overall size of the wavelet on the 
sky ^ MW (a) « 2 V2 a, whereas the second orthogonal size defines 
the size of the internal structure of the wavelet ^2 MW (^) ~ bn/ko 
(McEwen et al. 2005b). 



3 CROSS-CORRELATION IN WAVELET SPACE 

The effectiveness of using a spherical wavelet estimator to detect 
cross-correlations between the CMB and the nearby galaxy density 
distribution has been demonstrated by Vielva et al. (2006). We ex- 
tend the wavelet estimator here to account for directional wavelets 
that are not azimuthally symmetric. Firstly, we review the theo- 
retical cross-power spectrum of the CMB and galaxy density and 
give an expression for the spectrum based on the particular cosmo- 
logical model. We then define the wavelet coefficient covariance 
estimator used to test for cross-correlations. Using the theoretical 
cross-power spectrum previously described, we give a theoretical 
prediction for the wavelet coefficient covariance for a particular 
cosmological model (i.e. for a particular theoretical cross-power 
spectrum). In subsequent sections, we use this result to compare 
theoretical predictions of the wavelet covariance for a range of cos- 
mologies with measurements from the data, in order to place con- 
straints on the cosmological parameters that define the dark energy. 
Finally, we also examine in this section the effectiveness of various 
spherical wavelets for detecting cross-correlations. 

3.1 Theoretical cross-power spectrum 

For a particular cosmology, we consider the theoretical cross-power 
spectrum C^ T of the galaxy density map n{co) with the CMB tem- 




(a) Spherical Mexican hat wavelet (SMHW) 




(b) Spherical butterfly wavelet (SBW) 




(c) Spherical real Morlet wavelet (SMW) 



Figure 2. Spherical wavelets at scale a = b = 0.2. Wavelet maps are dis- 
played in the Mollweide projection, where the wavelets have been rotated 
down from the north pole for ease of observation. The SMW is plotted for 
wave vector k = (10, 0) r . 

perature anisotropy map t(a>), defined as the ensemble average of 
the product of the spherical harmonic coefficients of the two maps: 

( n €m ffm') = <W<W Cf T , (7) 

where 5^ is the Kronecker delta function. In defining the cross- 
correlation in this manner we implicitly assume that the galaxy den- 
sity and CMB random fields on the sphere are homogeneous and 
isotropic. Modelling the observation process, the observed cross- 
spectrum for the given cosmology is related to the theoretical one 
by 

C^ T,obs = Pe 2 b*f b] Cf , (8) 

where is the pixel window function for the pixelisation scheme 
adopted and b^ and b T £ are the beam window functions for the 
galaxy density 7 and CMB maps respectively. Pixel window and 
beam functions are represented here by Legendre coefficients 
(rather than spherical harmonic coefficients). 

The theoretical cross-power spectrum may be computed for a 
given cosmology by (e.g. Nolta et al. 2004) 

Cf = 12/rQ m // 2 J y A*(*) Ff(k) F](k) , (9) 

where Q m is the matter density, H is the Hubble parameter, 
Ag(k) = k 3 Ps(k)/27T 2 is the logarithmic matter power spectrum 
(Ps(k) is the matter power spectrum) and Ff(k) and Fj(k) are the 

7 There is no beaming for the NVSS galaxy density map subsequently used, 
hence = 1, W. 
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filter functions for the galaxy density distribution and CMB respec- 
tively, given by 



Ff(k) = b 



D(z) jdhfe)] 



and 



Fj(k) 



Mkrj(z)] . 



(10) 



(11) 



The integration required to compute Fj(£) is performed over z 
from zero to the epoch of recombination, whereas, in practice, 
the integration range for F^(k) is defined by the source redshift 
distribution function The function D(z) is the linear growth 
factor for the matter distribution (calculated from CMBFAST 8 
(Seljak & Zaldarriaga 1996) by computing the transfer function for 
different redshifts), g = (1 + z)D(z) is the linear growth suppression 
factor, je[krj(z)] is the spherical Bessel function and t/(z) is the con- 
formal look-back time. The bias factor b is assumed to be redshift 
independent. 

For the evolution of D(z) we consider both the standard model 
dominated by a cosmological constant and also alternative models 
dominated by a time varying dark energy with negative pressure, 
p = p (l + z) 3(1+vv) , and energy density that is spatially inhomoge- 
neous. These alternative models are parameterised by the equation- 
of-state parameter w, defined as the ratio of pressure to density. 
The standard inflationary model assumes a homogeneous field (the 
cosmological constant) with w = -1, however in general dark en- 
ergy models are characterised by values of w < 0. For instance, 
topological defects can be phenomenologically represented by an 
equation-of- state parameter -2/3 < w < -1/3 (e.g. Friedland et al. 
2003), quintessence models imply -1 < w < (Wetterich 1988; 
Caldwell et al. 1998) and phantom models have equation-of- state 
parameter w < -1 (these last models, however, violate the null 
dominant energy condition; see Carroll et al. (2003) for a detailed 
discussion). 

A convincing explanation for the origin and nature of the dark 
energy within the framework of particle physics is lacking, never- 
theless in general the models considered in the literature produce a 
w that varies with time. For most dark energy models, however, the 
equation-of- state changes slowly with time and a standard approxi- 
mation is that (at least during a given epoch) w can be considered as 
a constant equation-of- state parameter (Wang et al. 2000). Hence- 
forth, we consider w to be constant as a useful approach to extract 
fundamental properties of the dark energy. 



3.2 Wavelet covariance estimator 

The covariance of the wavelet coefficients is used as an estimator to 
detect any cross-correlation between the CMB and the galaxy den- 
sity distribution. A positive covariance indicates a positive cross- 
correlation between the data. The wavelet coefficient covariance 
estimator is defined as the sum over all points on the wavelet do- 
main sky, of the product of the wavelet coefficient maps: 

b, y) = y <# W ^ b > a ^ y"> b > a >P> 7) ' ( 12 ) 



where N a p is the number of samples in the wavelet domain sky, v aJ3 
is a weighting function and W™(a, b, a,p,y) and Wl(a,b,a,/3,y) 
are the wavelet coefficients of the galaxy density distribution and 
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CMB respectively. We choose weights to reflect the relative size 
of the pixels on the sky in the wavelet domain (equiangular sam- 
pling is used in the wavelet coefficient Euler angle domain), where 
v a p = \ sin/?. This weighting scheme ensures regions near the poles 
of the coordinate system do not have a greater influence on the es- 
timated covariance than regions near the equator. 

We may also average the covariance estimator over orienta- 
tions, so that we obtain an overall covariance measure for the given 
scales: 



X^(aM=jrJ]x™(a,b,y), 



(13) 



where N 7 is the number of samples in the wavelet domain ori- 
entational component. This measure is still sensitive to direc- 
tional structure when using a directional spherical wavelet, just 
as X^ T (a, b,y) and X^ T (a, b) are both sensitive to localised spatial 
structure in (a,p). 

A theoretical prediction of the wavelet covariance may be 
specified for a given cosmological model. Vielva et al. (2006) de- 
rive this for azimuthally symmetric spherical wavelets, however 
the extension to directional wavelets is non-trivial. We present 
the derivation of the theoretical wavelet covariance for directional 
wavelets in Appendix A, stating here the expression for the theo- 
retical covariance obtained: 



Xf{a^y): 



(14) 



where i// £m are the spherical harmonic coefficients of the wavelet. 
In practice at least one of the functions in (14) has a finite band 
limit so that negligible power is present in those coefficients above 
a certain £ max . All summations over £, here and subsequently, may 
therefore be truncated to £ max . 



3.3 Comparison of wavelet covariance estimators 

We compare the expected performance of various spherical 
wavelets for detecting the ISW effect in this section. Vielva et al. 
(2006) perform a similar analysis to compare the performance of 
the SMHW estimator to real and harmonic space estimators, show- 
ing the effectiveness of the SMHW estimator. Instead we focus here 
on comparing wavelets; thus we extend the analysis to directional 
wavelets. 

The expected signal-to-noise ratio (SNR) of the wavelet co- 
variance estimator of the CMB and the LSS density distribution 
may be used to compare the expected performance of various 
wavelets for detecting the ISW effect. The expected SNR for a 
particular scale is given by the ratio of the expected value of the 
wavelet covariance estimator and its dispersion: 



SNR^(<2,&) : 



&Xf(a,b) 



(15) 



where for a directional wavelet the variance of the wavelet covari- 
ance estimator is given by 

[AX™(a,b)f = ([x™(a,b)f)-(xf(a,b)) 2 

x |X |(^) fol | 2 ) ((C? T ) 2 + Q"cr) ,(16) 
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where Cj T and C^ N are the CMB and galaxy count power spectra 
respectively. For the case of azimuthally symmetric wavelets (16) 
reduces to the form given by Vielva et al. (2006). 9 

The expected SNR is computed for the spherical wavelets we 
consider for a range of dilations (the dilations considered are those 
that are subsequently used to attempt to detect the ISW effect, as 
defined in Table 1). The theoretical power spectra used in this ex- 
periment are computed from the cosmological concordance model 
parameters specified in Table 1 of Spergel et al. (2003) (£2 A = 0.71, 
Q m = 0.29, Q b = 0.047, H = 72, r = 0.166, n = 0.99). CMBFAST 
is used to simulate the CMB spectrum Cj T , the cross-power C^ T 
spectrum is simulated in accordance with the theory outlined in sec- 
tion 3.1, while the actual NVSS data is used to provide the LSS an- 
gular power spectrum C^ N . Ideal noise-free and full-sky conditions 
are assumed for this experiment. Beam and pixel windowing are not 
considered in the results presented, although when included these 
factors make only a very minimal difference to the numerical values 
obtained. The expected SNR computed for each wavelet for a range 
of scales is illustrated in Fig. 3. The relative expected performance 
of the SMW for detecting the ISW effect is low, thus we do not 
consider this wavelet any further. The SMHW and the SBW are of 
comparable expected performance. However, notice that the ellipti- 
cal SMHW (illustrated by the off-diagonal values of Fig. 3 (b)) and 
the SBW are slightly superior to the azimuthally symmetric SMHW 
(illustrated by the diagonal values of Fig. 3 (b)). It may appear 
counter intuitive initially that asymmetric wavelets may perform 
better at detecting the ISW effect when the data are isotropic fields. 
Although the statistics of the fields are globally isotropic, this does 
not preclude local oriented features in the data. Indeed, the distri- 
bution of the ellipticity of peaks in a isotropic Gaussian random 
field has been derived by Barreiro et al. (1997, 2001), illustrating 
that one would expect to see local rotationally invariant features. 
The maximum expected SNR achievable with both the SMHW and 
SBW is approximately 4, compared to an estimated best achievable 
SNR of approximately 7.5 for a perfect survey (Afshordi 2004). 



4 ANALYSIS PROCEDURE 

We describe in this section the data and an overview of the analysis 
procedure used to attempt to detect the ISW effect using directional 
spherical wavelets. 

4.1 WMAPdata 

To minimise the contribution of foregrounds and systematics to 
CMB anisotropy measurements the WMAP assembly contains a 
number of receivers that observe at a range of frequencies. The 
WMAP team and other independent groups have proposed vari- 
ous constructions of CMB maps from data measured by differ- 
ent receivers and bands in order to minimise foreground contri- 
butions (template based methods: Bennett et al. 2003b; template 
independent methods: Bennett etal. 2003b, Tegmarketal. 2003, 
Eriksen et al. 2004). In this work, we use the template based fore- 
ground removed maps (Bennett et al. 2003b) to construct the co- 
added map proposed by the WMAP team and used in their non- 
Gaussianity analysis (Komatsu et al. 2003). Following the data pro- 
cessing pipeline specified by Komatsu et al. (2003), the foreground 

9 Note that Vielva et al. (2006) use Legendre coefficients to represent the 
symmetric SMHW, rather than spherical harmonic coefficients. 
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(c) SBW for all dilations 

Figure 3. Expected SNR of the wavelet covariance estimator of CMB and 
radio source maps (cosmological parameters are chosen according to the 
concordance model). All spherical wavelets are considered for isotropic 
dilations in panel (a): SMHW (blue, dashed); SBW (red, dotted-dashed); 
SMW (magenta, dotted). Due to the relatively poor expected performance 
of the SMW, the SNR for all anisotropic dilations in only shown for the 
SMHW and SBW in panels (b) and (c) respectively. Contours at SNR val- 
ues of two, three and four are also plotted in panels (b) and (c). 



cleaned WMAP maps for which the CMB is the dominant signal 
(two Q-band maps at 40.7GHz, two V-band maps at 60.8GHz and 
four W-band maps at 93.5GHz) are combined to give the signal- 
to-noise ratio enhanced co-added map. The conservative KpO ex- 
clusion mask provided by the WMAP team is applied to remove 
remaining Galactic emission and bright point sources. The fore- 
ground removed maps used to construct the co-added map and the 
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(b) NVSS 

Figure 4. WMAP and NVSS maps after application of the joint mask. The 
maps are represented in the HEALPix format at a resolution of = 64 
(corresponding to a pixel size of -55'). 

KpO mask are available from the Legacy Archive for Microwave 
Background Data Analysis (LAMBDA). 10 The maps are provided 
in the HEALPix 11 (Gorski et al. 2005) format at a resolution of 
Mide = 512 (giving \2N^ Q 2 ~ 3 x 10 6 pixels on the sphere). 
Since the ISW signal we hope to detect is expected to be generated 
by structure on scales greater than 2° (Afshordi 2004), we down- 
sample the co-added map constructed to Af side = 64 (~ 5 x 10 4 
pixels). This corresponds to a pixel size of approximately 55' and 
so should be of sufficient resolution to detect the ISW signal. The 
co-added WMAP map used in the subsequent analysis is displayed 
in Fig. 4 (a) (with the joint KpO-NVSS mask applied, as described 
in section 4.2). 

4.2 NVSS data 

We use the NVSS radio source catalogue as a local tracer of the 
LSS. The catalogue covers approximately 80% of the sky and con- 
tains measurements of almost 2 x 10 6 point sources with a mini- 
mum flux density of ~2.5mJy. Although the distance of individ- 
ual sources is largely unknown, the sources are thought to be dis- 
tributed in the redshift range < z < 2 with a peak in the distri- 
bution at z ~ 0.8 (Boughn & Crittenden 2002). The ISW signal we 
hope to detect is expected to be produced at z ~ 0.4 with negligible 
contribution from z > 1.5 (Afshordi 2004). Furthermore, the corre- 
lations induced by the ISW effect are a relatively large scale phe- 
nomenon, hence any detection will be cosmic variance limited. The 
near full- sky coverage and source distribution of the NVSS data 
thus make it an ideal probe of the local matter distribution to use 
when searching for the ISW effect. Indeed, the NVSS catalogue has 

10 http://cmbdata.gsfc.nasa.gov/ 

11 http://healpix.jpl.nasa.gov/ 



been used to make positive detections of the ISW effect already us- 
ing the real space correlation function (Boughn & Crittenden 2002, 
2004; Nolta et al. 2004). 

As first noticed by Boughn & Crittenden (2002), the NVSS 
catalogue mean source density varies with declination. We correct 
this systematic by scaling each iso-latitude band to impose a con- 
stant mean density. The iso-latitude bands are defined by the con- 
tinuous root-mean- squared (RMS) noise regions of the survey (see 
Fig. 10 of Condon et al. 1998). Alternative correction strategies 
were also considered (e.g. the one proposed by Nolta et al. 2004), 
however we found cross-correlation detections to be insensitive to 
the particular correction procedure. 

We represent the corrected NVSS source distribution in the 
HEALPix format at the same resolution of the WMAP data 
Wide = 64). We consider only the sources in the NVSS catalogue 
above 2.5mJy, corresponding to an identification completeness of 
50% (Condon et al. 1998). No observations are made in the NVSS 
catalogue for equatorial declination 5 lower than -50° and the cov- 
erage in the range -50° < 6 < -37° is not sufficient. Hence, we 
consider only sources with an equatorial declination 5 ^ -37°. 
With these constraints we obtain a galaxy distribution map con- 
taining ~1.6 x 10 6 sources with an average number of 40.4 counts 
per pixel. Since the WMAP mask and NVSS coverage each exclude 
various sections of the sky, we construct a joint mask (KpO + 6 < 
-37°) that leaves only those common pixels that remain in both 
maps. This mask is applied to both the co-added WMAP and NVSS 
maps used in the subsequent analysis. The NVSS map with the joint 
map applied is displayed in Fig. 4 (b). 

4.3 Simulations 

Monte Carlo simulations are performed to construct significance 
measures for the wavelet covariance statistics used to detect the 
ISW effect. 1000 Gaussian simulations of the WMAP data are 
constructed from the power spectrum produced by CMBFAST us- 
ing the cosmological concordance model parameters specified in 
Table 1 of Spergel et al. (2003). For each realisation we simu- 
late the WMAP observing strategy and then construct a simulated 
co- added map. Measurements made by the Q-, V- and W-band 
receivers are simulated by convolving with realistic beams and 
adding anisotropic WMAP noise for each receiver. The procedure 
described in section 4.1 to construct the co-added map is then per- 
formed on the simulated data. The simulated co-added map is then 
down- sampled to A^ sid e = 64 and the joint mask is applied. 

4.4 Procedure 

The analysis procedure consists of computing the wavelet covari- 
ance estimator described in section 3.2 from the wavelet coeffi- 
cients of the co- added WMAP and NVSS data computed for a 
range of scales and, for the directional wavelets, a range of y orien- 
tations. We consider only those scales where the ISW signal is ex- 
pected to be significant (Afshordi 2004): the wavelet scales consid- 
ered and the corresponding effective size on the sky of the wavelets 
are shown in Table 1 . Each of the directional wavelets considered 
are rotationally invariant under integer azimuthal rotations of tt, 
thus the azimuthal rotation angle y effectively lies in the domain 
[0,7r). For directional wavelets we consider five evenly spaced y 
orientations in the domain [0,tt). Any deviation from zero in the 
wavelet covariance estimator for any particular scale or orientation 
is an indication of a correlation between the WMAP and NVSS 
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data and hence a possible detection of the ISW effect. An identical 
analysis is performed using the simulated co-added CMB maps in 
place of the WMAP data in order to construct significance mea- 
sures for any detections made. Finally, we use any detections of the 
ISW effect to constrain dark energy parameters. 

The application of the joint KpO-NVSS mask distorts wavelet 
coefficients corresponding to wavelets with support that overlaps 
the mask exclusion region. These contaminated wavelet coefficients 
must be removed from the wavelet covariance estimate. We con- 
struct an extended coefficient mask for each scale to remove all 
contaminated wavelet coefficients from the analysis. Binary mor- 
phological operations, commonly applied on the plane in image 
processing (Gonzalez 2003), are extended to the sphere and applied 
to construct the extended coefficient mask. A morphological dila- 
tion adds pixels to the boundaries of regions, whereas a morpho- 
logical erosion removes pixels from boundaries. Firstly, the orig- 
inal joint mask is 'opened' (a morphological dilation followed by 
an erosion) to remove point source regions whilst maintaining the 
size of the original central masked region. The central masked re- 
gion is then extended by half the effective size of the wavelet (at the 
particular scale) by performing a morphological erosion. The orig- 
inal point source mask regions are then replaced by applying the 
original mask (point source regions are removed initially to ensure 
that they are not expanded when performing the final erosion). The 
result of this procedure is to extend the central region of the joint 
mask by half the effective size of the wavelet, while maintaining the 
point source regions. This ensures that all coefficients correspond- 
ing to wavelets that overlap with the masked Galactic plane are 
excluded from the analysis. However, wavelets that overlap with 
masked point sources are not removed. Retaining these wavelet co- 
efficients induces minimal distortion due to the large support of the 
wavelets relative to the size of the masked point sources. This is 
a less conservative procedure for constructing extended coefficient 
masks than used in our non-Gaussianity analysis (McEwen et al. 
2005b), however the low resolution of the data and the large initial 
joint mask reduce drastically the number of pixels available ini- 
tially. If we construct more conservative extended masks we find 
that we quickly and significantly reduce the number of pixels con- 
sidered and thus the efficiency of the analysis. To ensure any cor- 
relations detected are not due to contamination that is not excluded 
by the extended mask we have also performed a number of tests 
using masks with larger exclusion regions. We find that the cor- 
relation signals that we subsequently detect (see section 5) remain 
when using these larger exclusion masks. Hence our choice of mask 
appears appropriate as it is more efficient than more conservative 
masks but still removes the contaminated wavelet coefficients. 



5 RESULTS AND DISCUSSION 

In this section we describe the results obtained when performing 
the directional spherical wavelet analysis procedure described pre- 
viously. We describe the detections of the ISW effect that we make, 
before using the detections to constrain dark energy parameters. 

5.1 Detection of the ISW effect 

Firstly, we describe the detection of positive correlations between 
the WMAP and NVSS data using various wavelets. To examine 
the source of the correlation detected we perform some prelimi- 
nary tests to ascertain whether perhaps foregrounds or systematics 
are responsible. We then localise regions on the sky that contribute 



most strongly to the correlation signal and examine these in more 
detail. 



5.1.1 Detections 

The wavelet covariance of the WMAP-NVSS data are shown in 
Fig. 5 for each of the wavelets considered. The covariance values 
are shown in units of N^, the number of standard deviations that 
the data measurement differs from the mean of 1000 Monte Carlo 
simulations for each covariance statistic. On examining the distri- 
bution of the wavelet covariance statistics of the simulations, the 
covariance statistics themselves appear to approximately Gaussian 
distributed. This implies that the approximate significance of any 
detection of a non-zero covariance can be inferred directly from the 
No- level. 12 In any case, 68%, 95% and 99% significance regions are 
plotted in Fig. 5 also (these significance levels are computed from 
the simulations directly but correspond approximately to l.OOcr, 
1.96cr and 2.51 cr respectively). The left panel in Fig. 5 shows the 
results for the azimuthally symmetric SMHW (with isotropic dila- 
tions), whereas the other panels correspond to directional wavelets. 
These latter panels exhibit a jagged structure due to the ordering 
of the wavelet scale index the covariance is plotted against. The 
wavelet scale index is ordered by sequentially concatenating rows 
from the 2- dimensional (a, b) dilation space. It would be preferable 
to view the data as a surface in (a, b) space, however it would not be 
possible to view the data and the six significance level surfaces at 
once. The scale index j used in Fig. 5 and subsequently may be re- 
lated to the scale index defined in Table 1 for a and b by |_(/ - 1)/7J 
and mod(j* - 1,7) + 1 respectively, where |_*J is the largest integer 
less than x and mod(-,/i) is the modulo n function. For all of the 
wavelets considered, the wavelet covariance of the data lies outside 
of the 99% significance level on certain scales. 13 

In Fig. 6 we show the No- surfaces for each wavelet in (a, b) 
space. The maximum N a of all detections made with the directional 
SMHW and SBW is 3.9 for both wavelets (the maximum detec- 
tion made with the SBW is slightly greater than that made with 
the SMHW, however to one decimal place they are equivalent). 
Using the symmetric SMHW (i.e. using only isotropic dilations) 
however, one obtains a maximum No- of 3.6. It should be noted that 
the use of the maximum No- to characterise the detection is based 
on an a posteriori scale selection: the significance of the detection 
is characterised for a specific scale, rather than the range of scales 
examined. It is interesting to note that for both wavelets the No- sur- 
faces displayed in Fig. 6 are perfectly consistent with the expected 
SNR displayed in Fig. 3: the most significant detections are made 
on scales with high expected SNR. Note that the No- surface for the 
SMHW (Fig. 6 (a)) is not perfectly symmetric about the line a = b. 
Although the wavelet has identical shape for dilations (a, b) and 
(b,a), the second wavelet is rotated azimuthally by tt/2 relative to 
the first wavelet. Since the starting position differs between the two 

12 Ideally one would compute the significance of a detection directly from 
the Monte Carlo simulations. Indeed we do this to construct significance 
levels in Fig. 5, however for 1000 simulations the maximum significance 
one can measure is 99.9%. For many cases this does not give a fine enough 
resolution, thus we use the No- level to compare the significance of detec- 
tions made on different scales and with different wavelets. For computa- 
tional reasons we cannot extend easily the number of simulations far above 
1000. 

13 Although the SMW was shown to be poor at detecting cross-correlations 
in the data we ran the analysis regardless. As expected, the analysis results 
showed that the wavelet is indeed ineffective at detecting any correlation. 
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Table 1. Wavelet scales considered in attempting to detect the ISW effect. The effective sizes on the sky of 
the various wavelets for each wavelet scale are also shown. Note that for anisotropic dilations the wavelets 
have two effective sizes on the sky defined in orthogonal directions (see text). 
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(a) Symmetric SMHW (b) Elliptical SMHW (c) SBW 

Figure 5. Wavelet covariance statistics in units of No- for the WMAP and NVSS data against wavelet scale index. Significance levels obtained from 1000 
Monte Carlo simulations are shown by the shaded regions for 68% (yellow/light-grey), 95% (magenta/grey) and 99% (red/dark-grey) levels. See text for a 
discussion of the ordering of the dilation index. 



cases the directional analysis of each probes slightly different di- 
rections, resulting in an No- surface that is not perfectly symmetric 
about the diagonal. 

The symmetric SMHW analysis has already been performed 
by Vielva et al. (2006). We have repeated the analysis here as a 
consistency check for our analysis and also to confirm the previous 
work. The wavelet covariance curve shown in Fig. 5 (a) is similar 
to the curve obtained by Vielva et al. (2006), however we make a 
slightly more significant detection at 3.6<x compared to the detec- 
tion made by Vielva et al. (2006) at 3.3cr. After performing numer- 
ous tests with different weighting and masking schemes we believe 
the difference in significance with Vielva et al. (2006) is due to the 
different mask construction techniques adopted, the different pix- 
elisation of the wavelet domain and the slightly different definition 
of the SMHW used (see Martinez-Gonzalez et al. (2002) for the 
definition used by Vielva et al. (2006)). 

In summary, we have detected a positive correlation between 
the WMAP and NVSS data, as indicated by a positive wavelet 
covariance, at the 3.9<r level on wavelet scales about (a,b) = 
(100', 300'). The sign of the correlation detected (positive) and the 
scale the detections are made at is consistent with an ISW sig- 
nal. We next perform some preliminary tests to determine whether 
the correlation we detect is indeed produced by the ISW effect or 
whether other factors are responsible. 



5.1.2 Foregrounds and systematics 

To test whether foregrounds or WMAP systematics are responsi- 
ble for the correlation signal we examine the wavelet covariance 
obtained using both the separate WMAP bands and also combina- 
tions of difference maps constructed from the individual bands. We 
have two Q-band maps at 40.7GHz, two V-band maps at 60.8GHz 
and four W-band maps at 93.5GHz to construct test maps from (we 
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Figure 6. Wavelet covariance No- surfaces. Contours are shown for levels 
of two and three No-. 
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subsequently use the notation Ql, Q2, etc. to denote the signal mea- 
sured by each receiver and the single band letter to denote the sum 
of all maps measured for the given band, e.g. Q=Q1+Q2). Firstly, 
the wavelet covariance of the NVSS map with each of the Q-, V- 
and W-band maps is computed and is displayed in Fig. 7. For each 
case the wavelet covariance signal is essentially identical to the sig- 
nal observed when using the co-added WMAP map (which is also 
plotted in Fig. 7 for comparison). Secondly, we construct the dif- 
ference maps Q1-Q2, V1-V2 and W1-W2+W3-W4 for each band 
from the signals measured by different receivers. The CMB signal 
and foreground contributions are completely removed from these 
noisy difference maps. In addition, we examine the difference map 
W-V-Q which is free of CMB but has a clear foreground contribu- 
tion. The wavelet covariance of the NVSS map with each of these 
difference maps is computed and is displayed in Fig. 8. For all of 
these difference maps the correlation previously detected in the co- 
added WMAP is eliminated. 

Any correlation between the WMAP and NVSS maps due to 
unremoved foreground emission in the WMAP data is expected to 
be frequency dependent, reflecting the emission law of the point 
source population (e.g. Toffolatti et al. 1998). However we observe 
identical wavelet covariance signals in each band (Fig. 7) and do 
not detect any frequency dependence. Furthermore, the W-V-Q dif- 
ference map has minimal CMB contribution but has a clear fore- 
ground contribution. We do not detect any correlation using this 
map (Fig. 8). These findings suggest that the correlation signal that 
we detect in the co-added WMAP map is not due to unremoved 
foregrounds in the WMAP data (we show also in section 5.1.3 that 
the detection is not due to a few localised regions, giving further 
evidence to support the claim that unremoved foregrounds are not 
responsible). 

The wavelet covariance signal that we detect in the co- added 
WMAP map is present in all of the individual WMAP bands 
(Fig. 7). Furthermore, the covariance signal is eliminated in each 
of the individual band difference maps that contain no CMB or 
foreground contributions (Fig. 8). Since all of the WMAP receivers 
produce an identical covariance signal that does not appear to be 
due to the noise artifacts of each receiver, we may conclude that it 
is unlikely that systematics are responsible for the detection made. 

These preliminary tests show that it is unlikely that unremoved 
foreground emission or systematic effects in the WMAP data are 
responsible for the detection of a positive correlation between the 
WMAP and NVSS data. It would appear that the detection is due to 
solely the CMB contribution of the WMAP data. This gives further 
support to the claim that the correlation we observe is produced by 
the ISW effect. 
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Figure 7. Wavelet covariance statistics in units of No- computed for the 
NVSS map with individual WMAP band maps. 




— I — W-V-Q 
Q1-Q2 
V1-V2 

W1-W2+W3-W4 



5 10 15 20 25 30 35 40 45 

Dilation index 

(a) SMHW 



5.1.3 Localised regions 

A wavelet analysis inherently allows the spatial localisation of in- 
teresting signal characteristics. The regions that contribute most 
strongly to the wavelet covariance detected may therefore be lo- 
calised, not only in scale, but also in position and orientation on the 
sky. 

The wavelet covariance statistic is essentially the weighted 
mean of the wavelet domain product map constructed from the 
product of the WMAP and NVSS wavelet coefficients. The regions 
that contribute most strongly to the covariance may thus be lo- 
calised by thresholding the wavelet domain product map. In Fig. 9 
we illustrate the localised regions detected once the wavelet prod- 
uct map is thresholded so that only those values that lie above 3cr 
remain. For illustrational purposes we show for each wavelet only 
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Figure 8. Wavelet covariance statistics in units of No- computed for the 
NVSS map with WMAP band difference maps. 
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the localised maps corresponding to the scale and first orientation 
of the most significant covariance detection. In the subsequent anal- 
ysis all orientations are considered, thus Fig. 9 is illustrative rather 
than representing all of the localised regions that we detect. No- 
tice that the localised regions found are relatively evenly distributed 
over the sky (remember that we cannot search within the extended 
joint WMAP-NVSS mask, which becomes progressively larger on 
larger wavelet scales). The localised maps found with the various 
wavelets also exhibit a fair degree of similarity. 

A interesting next step is to determine whether these localised 
regions are the sole source of the wavelet covariance signal de- 
tected. To examine this hypothesis we remove the localised re- 
gions detected on all scales and orientations (not just the ones 
corresponding to the maximum detections that are displayed in 
Fig. 9) and repeat the analysis. The wavelet covariance signals mea- 
sured with and without the localised regions removed are shown in 
Fig. 10. Once the localised regions are removed the covariance sig- 
nal is reduced in significance, as one would expect since those re- 
gions that contribute most strongly are removed, but the covariance 
signal is not eliminated. Indeed, the signal remains on certain scales 
at almost the 95% significance level. We conclude that the localised 
regions that we detect are not the sole source of correlation between 
the WMAP and NVSS data. Our findings are consistent intuitively 
with the predicitions of the ISW effect, namely that we would ex- 
pect any observed correlation to be due to weak correlations over 
the entire sky rather than a few localised regions. Furthermore, the 
aggregate nature of the signal detected provides further evidence to 
suggest that foreground contamination in the data is not responsi- 
ble. Boughn & Crittenden (2005) make similar conclusions that the 
correlation signal they detect is due to aggregate correlations over 
the entire sky rather than a few localised regions. 

As a final test of the integrity of the data we examine in more 
detail a few of the most significant localised regions that we detect. 
We examine 18 of the localised regions illustrated in Fig. 9, se- 
lected from the regions that contain the greatest wavelet coefficient 
product map values, in both the NVSS and Bonn 1420MHz (Reich 
1982; Reich & Reich 1986) radio surveys (the latter data set affords 
better visualisation of the LSS and enables one to highlight regions 
in the NVSS data for closer examination). The regions selected are 
defined in Table 2 and Fig. 1 1 . We do not notice any difference in 
the NVSS data between our localised regions and regions selected 
at random. The localised regions we detect do not, in general, cor- 
respond to regions with particularly bright point sources that could 
have caused contamination of the data. 

All of the tests that we have performed indicate that the cor- 
relation detected is indeed due to the ISW effect and not due to 
systematics or foreground contamination. Next we use our positive 
detection of the ISW effect to place constraints on dark energy. 

5.2 Constraints on £l A and w 

In a flat universe the ISW effect is present only in the presence 
of dark energy. We may therefore use our detection of the ISW 
effect to constrain dark energy parameters by comparing the theo- 
retically predicted wavelet covariance signal for different cosmo- 
logical models with that measured from the data. In particular, 
we constrain the vacuum energy density £l A and the equation- 
of-state parameter w (assuming the equation-of- state parameter 
does not evolve with redshift within the epoch of interest). We 
probe the parameter ranges < £l A < 0.95 and -2 < w < 0. 
For the other cosmological parameters we use the concordance 
model values (Spergel et al. 2003). We use a bias parameter value 



(a) Symmetric SMHW; dilation (a,b) = (250', 250') 
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Figure 9. Localised WMAP-NVSS wavelet coefficient product maps 
thresholded at 3cr. The regions remaining show those areas that contribute 
most strongly to the positive wavelet covariance signal detected. Localised 
product maps are only shown for the first orientation for the scale of the 
most significant correlation detection made with each wavelet. 




Figure 11. Approximate localised regions flagged for closer examination 
(see Table 2 for more details). 



of b = 1.6. This is the best choice consistent with the concor- 
dance model that is within the range favoured by the NVSS data 
(Boughn & Crittenden 2002, 2004) when adopting the RLF1 model 
proposed by Dunlop & Peacock (1990). In any case, Vielva et al. 
(2006) find dark energy constraints to be insensitive to changes 
in the bias within the range 1.4 < b < 1.8. We also use the 
Dunlop & Peacock (1990) RLF1 model to describe the ^ distri- 
bution since it fits the NVSS galaxy autocorrelation function well, 
as previously shown by Boughn & Crittenden (2002, 2004) and 
Nolta et al. (2004). 

To compare the wavelet covariance of the data with the theo- 
retical predictions made by different cosmological models we com- 
pute thex 2 f° r parameters = (£2 A , w), defined by 

X H@) = [x™(&) - Xff C" 1 [X* T (0) - Xf] , (17) 
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(a) Symmetric SMHW 



(b) Elliptical SMHW 



(c) SBW 



Figure 10. Wavelet covariance statistics computed for the WMAP and NVSS data with (green/light-grey dots) and without (blue/dark-grey pluses) the localised 
regions that contribute most strongly to the covariance removed. The covariance signal is reduced in significance when the localised regions are removed, as 
one would expect, but it is certainly not eliminated suggesting that the localised regions are not the sole source of the wavelet covariance. 



Table 2. Approximate localised regions flagged for closer examination (see 
Fig. 11 for a map of these regions). The wavelets that flag each region are 
listed also (the same regions are flagged by both the symmetric and elliptical 
SMHW). Although these regions are specified precisely they represent only 
the approximate position of the regions localised by the wavelet analysis, 
which is often of the order of a number of degrees. 



Region 


Location 


Flagged by wavelet 




Longitude 


Latitude 


SMHW 


SBW 


1 


75° 


57° 


y 


X 


2 


75° 


53° 


X 


y 


3 


323° 


56° 


y 


X 


4 


321° 


62° 


X 


y 


5 


267° 


50° 


y 


X 


6 


268° 


45° 


X 


y 


7 


213° 


40° 


y 


X 


8 


223° 


30° 


X 


y 


9 


160° 


26° 


y 


y 


10 


94° 


-28° 


y 


X 


11 


81° 


-34° 


y 


X 


12 


118° 


-42° 


X 


y 


13 


20° 


-48° 


y 


X 


14 


34° 


-31° 


X 


y 


15 


230° 


-68° 


y 


X 


16 


204° 


-56° 


y 


y 


17 


186° 


-54° 


y 


y 


18 


218° 


-33° 


y 


X 



where Xj^ T (0) is the predicted covariance vector for the given cos- 
mological model and X^ T is the wavelet covariance vector of the 
data. These vectors are constructed by concatenating the wavelet 
covariance values computed for the considered dilations and orien- 
tations into vectors. The matrix C is a similarly ordered covari- 
ance matrix of the wavelet covariance statistics, computed from 
the 1000 Monte Carlo simulations. Making a Gaussian approxi- 
mation for the likelihood, we may write the likelihood function 
as X(0) °c exp [-^ 2 (0)/2j. The likelihood surfaces for the pa- 
rameter ranges discussed for both of the SMHW and the SBW 
are illustrated in Fig. 12. Marginalised likelihood distributions are 
also shown. Parameter estimates are made from the mean of the 
marginalised distributions, with 68%, 95% and 99% confidence 
regions constructed to ensure the required probability is met in 
the tails of the distribution. The parameter estimates and errors 
made with the various wavelets are as follows: Q A = 0.63^ }^, 



w = -0.77^ using the SMHW; Q A = 0.52+g;|g, w = -0.73^ 
using the SBW. Within error bounds these parameter estimates are 
consistent with one another and also with estimates made from 
numerous other analysis procedures and data sets. The errors we 
obtain on the parameter estimates are relatively large. Although 
wavelets perform very well when attempting to detect the ISW ef- 
fect since one may probe different scales, positions and orienta- 
tions, once all information is incorporated to compute the likeli- 
hood surface the performance of a wavelet analysis is comparable 
to other linear techniques, as expected 14 (see Vielva et al. (2006) 
for a detailed discussion and comparison of real, harmonic and (az- 
imuthally symmetric) wavelet space techniques for detecting cor- 
relations due to the ISW effect). 

5.3 Constraints on Q A 

It is also interesting to use our detection of the ISW effect to pro- 
vide an estimate of the overall evidence for the existence of dark 
energy when assuming a pure cosmological constant. We compute 
the likelihood distribution of just the vacuum energy density £l A in 
the presence of a cosmological constant by repeating the procedure 
described in section 5.2 for = £l A and w = -1. This essentially 
reduces to taking a slice through the 2-dimensional likelihood sur- 
faces shown in Fig. 12 at w = -1. The corresponding likelihood 
distributions for each wavelet are shown in Fig. 13. From the mean 
of the distributions we obtain the following vacuum energy den- 
sity estimates: Q A = 0.70^ using the SMHW; Q A = 0.57^;}* 
using the SBW. Within error bounds these parameter estimates are 
again consistent with each other and with estimates made using 
other techniques and data sets. 

We show also in Fig. 13 the cumulative probability P(£l A > x). 
One may read directly off this plot the evidence for the existence of 
dark energy above a certain level. For instance, using the SMHW 
we may state that Q A > 0.1 at the 99.96% level. Using the SBW 
we may state that Q A > 0.1 at the 99.71% level. The cummulative 
probability function falls off faster for the SBW, indicating that the 
SMHW happens to be more effective at ruling out a low vacuum 
energy density (i.e. the SMHW predicts a higher £l A ). In any case, 
we have very strong evidence for dark energy in the case of a pure 
cosmological constant. 



The results of different linear analyses, however, will not be exactly iden- 
tical due to differences in masking. 



Detection of the ISW effect and corresponding dark energy constraints 13 




(a) Full likelihood surface obtained using SMHW 



(b) Full likelihood surface obtained using SBW 




(c) Marginalised distribution for £l A obtained using SMHW (d) Marginalised distribution for £l A obtained using SBW 




-0.8 -0.6 -0.4 -0.2 



(e) Marginalised distribution for w obtained using SMHW 



(f) Marginalised distribution for w obtained using SBW 



Figure 12. Likelihood surfaces for parameters @ = (Q A ,w). The full likelihood surfaces are shown in the top row of panels, with 68%, 95% and 99% 
confidence contours also shown. Marginalised distributions for each parameter are shown in the remaining panels, with 68% (yellow/light-grey), 95% (light- 
blue/grey) and 99% (dark-blue/dark-grey) confidence regions also shown. The parameter estimates made from the mean of the marginalised distribution is 
shown by the triangle and dashed line. 



6 CONCLUSIONS 

We have performed a directional spherical wavelet analysis to 
search for correlations between the WMAP and NVSS data aris- 
ing from the ISW effect. Wavelets are an ideal tool for searching 
for the ISW effect due to the localised nature of the effect and the 
localisation afforded by a wavelet analysis. Furthermore, we use di- 
rectional spherical wavelets that allow one to probe not only scale 
and spatially localised structure but also orientated structure in the 
data. 

Using this technique we have detected a significant correla- 
tion between the WMAP and NVSS data. The detection is made 
using both the directional SMHW and the SBW, on approximately 
the same wavelet scale of (a, b) = (100', 300'), corresponding to 
an overall effective size on the sky of approximately 10°. The scale 
and positive sign of the detected correlation is consistent with the 
ISW effect. To test this hypothesis further we examined other pos- 
sible causes of the signal and after performing some preliminary 
tests ruled out both foreground contamination and systematics in 
the WMAP data. A wavelet analysis inherently provides localised 
information, thus we were able to localise those regions on the sky 
that contribute most strongly to the overall correlation. On remov- 
ing these regions from the analysis the correlation originally de- 
tected is reduced in significance, as one would expect, but it is 
not eliminated, remaining at the 95% significance level in some 
cases. It therefore appears that although the localised regions that 
we detected are a strong source of correlation, they are not the sole 
source. This is again consistent with predictions made using the 
ISW effect, where one would expect to observe weak correlations 
over the entire sky rather than just a few localised regions. Finally, 



we examined 18 of the most significant localised regions in more 
detail in high-resolution maps of the NVSS and Bonn 1420MHz ra- 
dio surveys. These regions are not atypical and do not, in general, 
correspond to regions with particularly bright point sources that 
could have caused contamination of the data. All of these findings 
suggest that the correlation that we have detected is due to the ISW 
effect. In our (nearly) flat universe the ISW effect is present only 
in the presence of dark energy, hence our detection of the effect is 
direct evidence for dark energy. 

We have used our detection of the ISW effect to constrain dark 
energy parameters by comparing theoretical predictions made by 
different cosmological models with the data. Constraining both the 
vacuum energy density £l A and the equation-of- state parameter w, 
we obtained the following parameter estimates: £l A = 0.63^ ^, 
w = -0.77^-35 using the SMHW; Q A = 0.52+°;*°, w = -0.73+°;^ 
using the SBW. We have therefore obtained independent estimates 
of the parameters that, within error bounds, are consistent with one 
another and also with estimates made from many other analysis 
procedures and data sets. The errors we obtain on the parameter 
estimates are relatively large. Although wavelets perform very well 
when attempting to detect the ISW effect since one may probe par- 
ticular scales, positions and orientations, once all information is 
incorporated when computing the parameter estimates the perfor- 
mance of the wavelet analysis is comparable to other linear tech- 
niques, as expected. We also constrain Q A for the case of a pure 
cosmological constant, i.e. w = -1. For this case we obtain the 
following parameter estimates: £l A = 0.70+° using the SMHW; 
£l A = 0.57+q jg using the SBW. Again the estimates are consis- 
tent with one another and with those made elsewhere. For this case 
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(a) Likelihood distribution for Qa obtained using SMHW 




(b) Likelihood distribution for obtained using SBW 




(c) Cumulative probability functions for P(£l\ > x) 

Figure 13. Likelihood distributions and cumulative probability functions 
for Qa when w = -1. Confidence regions at 68% (yellow/light-grey), 95% 
(light-blue/grey) and 99% (dark-blue/dark-grey) are also shown on the like- 
lihood distributions, with the parameter estimates made from the mean of 
the distribution shown by the triangle and dashed line. The cummulative 
probability functions show the probability P(Q\ > x) for the SMHW (solid 
magenta curve) and SBW (dashed blue curve). 



we also examined the evidence for a non-zero cosmological con- 
stant, ruling out £l A < 0.1 at greater than the 99% level using each 
wavelet. 

In this work we have extended to directional wavelets the 
spherical wavelet analysis technique used to probe CMB data and 
LSS tracers for cross-correlations. The effectiveness of the wavelet 
analysis technique is illustrated herein by the highly significant 
detection of the ISW effect that we are able to make using the 
WMAP and NVSS data. In addition, the spherical wavelet covari- 
ance approach to detecting correlations can be applied to other 
applications, such as searching for the Sunyaev-Zel'dovich effect 
(Hansen et al. 2005) or looking for common structures between the 
CMB and foregrounds (Liu & Zhang 2006). It would also be inter- 
esting to repeat our analysis using different tracers of the LSS and 
also the next release of the WMAP data to see if the ISW signal 
that we have detected remains. 



ACKNOWLEDGMENTS 

We thank Jacques Delabrouille and Daniel Mortlock for insightful 
discussions and comments on the manuscript and Dave Green for 
help accessing the NVSS and Bonn 1420MHz radio survey cata- 
logues. JDM would like to thank all at Universidad de Cantabria 
and College de France for the warm hospitality he received 
throughout his visit to each. JDM also thanks the Association of 
Commonwealth Universities and the Cambridge Commonwealth 
Trust for the support of a Commonwealth (Cambridge) Scholar- 
ship, King's College Cambridge for a Ferris Travel Grant and the 
Cambridge Philosophical Society for a travel grant. PV and EMG 
acknowledge financial support from the Spanish MEC project 
ESP2004-07067-C03-01. Some of the results in this paper have 
been derived using the HEALPix package (Gorski et al. 2005). We 
acknowledge the use of the Legacy Archive for Microwave Back- 
ground Data Analysis (LAMBDA). Support for LAMBDA is pro- 
vided by the NASA Office of Space Science. 



REFERENCES 

Afshordi N., 2004, Phys. Rev. D., 70, 83536 
Afshordi N., Loh Y. -S., Strauss M. A., 2004, Phys. Rev. D., 69, 
83524 

Antoine J. -P., Murenzi R., Vandergheynst P., Ali S. T., 2004, Two- 
dimensional wavelets and their relatives, Cambridge University 
Press, Cambridge 

Antoine J. -P. and Vandergheynst P., 1999, Applied and Computa- 
tional Harmonic Analysis, 7, 1 

Antoine J. -P. and Vandergheynst P., 1998, J. of Math. Phys., 39, 
8, 3987 

Antoine J. -P., Demanet L., Jacques L., 2002, Applied and Com- 
putational Harmonic Analysis, 13, 3, 177 

Barreiro R. B., Hobson M. P., Lasenby A. N., Banday A. J., 
Gorski K. M. and Hinshaw G., 2000, Mon. Not. Roy. Astron. 
Soc.,318, 475 

Barreiro R. B., Martmez-Gonzalez E., Sanz J. L., 2001, Mon. Not. 

Roy. Astron. Soc, 322, 411 
Barreiro R. B., Martmez-Gonzalez E., Sanz J. L., Cayon L., Silk 

J., 1997, ApJ, 478, 1 
Bennett C. L. et al. (WMAP team), 2003a, ApJS, 148, 1 
Bennett C. L. et al. (WMAP team), 2003b, ApJS, 148, 97 
Bennett C. L., Banday A. J., Gorski K. M., Hinshaw G., Jackson 

P., Keegstra, P., Kogut A., Smoot G. F, Wilkinson D. T., Wright 

E. L., 1996, ApJ, 464, 1 
Bogdanova I., Vandergheynst P., Antoine J. -P., Jacques L., Morvi- 

done M., 2005, Applied and Computational Harmonic Analysis, 

19, 2, 223 
Boldt E., 1987, Phys. Rep., 146, 215 

Boughn S. P., Crittenden R. G., 2005, Mon. Not. Roy. Astron. 

Soc, 360, 1013 
Boughn S. P., Crittenden R. G., 2004, Nature, 427, 45 
Boughn S. P., Crittenden R. G., 2002, Phys. Rev. Lett., 88, 21302 
Caldwell R. R., Dave R., Steinhardt P. J., 1998, Phys. Rev. Lett., 

80, 1582 

Carroll S. M., Hoffman M., Trodden M., 2003, Phys. Rev. D., 68, 
23509 

Cayon L., Jin J., Treaster A., 2005, Mon. Not. Roy. Astron. Soc, 
362, 826 

Cayon L., Martmez-Gonzalez E., Argueso F, Banday A. J., 
Gorski K. M., 2003, Mon. Not. Roy. Astron. Soc, 339, 1189 



Detection of the ISW effect and corresponding dark energy constraints 15 



Cayon L., Sanz J. L., Martmez-Gonzalez E., Banday A. J., 
Argueso R, Gallegos J. E., Gorski K. M., Hinshaw G., 2001, 
Mon. Not. Roy. Astron. Soc, 326, 1243 

Condon J. J. et al., 1998, ApJ, 115, 1693 

Corasaniti P. -S., Giannantonio T., Melchiorri A., 2005, 

Phys. Rev. D., 71, 123521 
Crittenden R. G., Turok N., 1996, Phys. Rev. Lett., 76, 575 
Cruz M., Tucci M., Martinez-Gonzalez E., Vielva P., 2006, Mon. 

Not. Roy. Astron. Soc., 369, 57 
Cruz M., Martmez-Gonzalez E., Vielva P., Cayon L., 2005, Mon. 

Not. Roy. Astron. Soc., 356, 29 
Dunlop S. J., Peacock J. A., 1990, Mon. Not. Roy. Astron. Soc, 

247, 19 

Eriksen H. K., Banday A. J., Gorski K. M., Lilje P. B., 2004, ApJ, 
612, 633 

Fosalba P., Gaztanaga E., 2004., Mon. Not. Roy. Astron. Soc, 
350, 37 

Fosalba P. Gaztanaga E., Castander E, 2003, ApJ, 597, 89 
Friedland A., Murayama H., Perelstein M., 2003, Phys. Rev. D., 
67, 43519 

Gonzalez R. L., 2001, Digital Image Processing, Second Edition, 
Prentice Hall 

Gorski K. M., Hivon E., Banday A. J., Wandelt B. D., Hansen F. 

K., Reinecke M., Bartelmann M., 2005, ApJ, 622, 759 
Hansen F. K., Branchini E., Mazzotta P., Cabella P., Dolag K., 

2005, Mon. Not. Roy. Astron. Soc, 361, 753 
Hu W., Scranton R., 2004, Phys. Rev. D., 70, 123002 
Jarrett T. H. et al., 2000, ApJ, 1 19, 2498 
Komatsu E. et al., 2003, ApJS, 148, 119 
Liu X., Zhang S. N., 2006, ApJL, 636, 1 
Liu X., Zhang S. N., 2005, ApJ, 633, 542 

Maddox S. J., Efstathiou G., Sutherland W. J., Loveday J., 1990, 

Mon. Not. Roy. Astron. Soc, 242, 43 
McEwen J. D., Hobson M. P., Mortlock D. J., Lasenby A. N., 

2005a, IEEE Trans. Sig. Proc, in press (astro-ph/0506308) 
McEwen J. D., Hobson M. P., Lasenby A. N., Mortlock D. J., 

2005b, Mon. Not. Roy. Astron. Soc, 359, 1583 
Martmez-Gonzalez E., Gallegos J. E., Argueso R, Cayon L., Sanz 

J. L., 2002, Mon. Not. Roy. Astron. Soc, 336, 22 
Mukherjee P., Wang Y, 2004, ApJ, 613, 51 
Nolta et al. (WMAP team), ApJ, 608, 10 

Padmanabhan N., Hirata C. M., Seljak U., Schlegel D. J., 
Brinkmann J., Schneider D. P., 2004, submitted to Phys. Rev. D., 
72, 043525 

Peebles P. J. E., Ratra B., 2003, Rev. Mon. Phys. 75, 599 
Perlmitter S. et al. (The Supernova Cosmology Project), 1999, 
ApJ, 517, 565 

Pogosian L., Corasaniti P. S., Stephan-Otto C, Crittenden R., 

Nichol R., 2005, Phys. Rev. D., 72, 103519 
Pogosian L., 2004, JCAP, 504, 15 
Reich W., 1982, A&AS, 48, 219 
Reich P., Reich W., 1986, A&AS, 63, 205 
Riess A. G. et al., 1998, Astron. J., 116, 1009 
Sachs R. K., Wolfe A. M., 1967, ApJ, 147, 73 
Scranton R. et al. (SDSS collaboration), 2003, submitted to 

Phys. Rev. Lett, (astro-ph/0307335) 
Seljak U., Zaldarriaga M., 1996, ApJ, 469, 437 
Spergel D. N. et al. (WMAP team), 2003, ApJS, 148, 175 
Tegmark M., de Oliveira-Costa A., Hamilton A. J. S., 2003, 

Phys. Rev. D., 68, 12 
Toffolatti L., Argueso Gomez E, de Zotti G., Mazzei P., Frances- 



chini, A., Danese, L., Burigana, C, 1998, Mon. Not. Roy. As- 
tron. Soc, 297, 117 
Tosic I., Frossard P., Vandergheynst P., 2005, Proc. IEEE DCC 
Vielva P., Martmez-Gonzalez E., Tucci M., 2006, Mon. Not. Roy. 

Astron. Soc, 365, 891 
Vielva P., Martmez-Gonzalez E., Barreiro R. B., Sanz J. L., Cayon 

L., 2004, ApJ, 609, 22 
Wandelt B. D., Gorski K. M., 2001, Phys. Rev. D., 63, 123002, 1 
Wang L., Caldwell R. R., Ostriker J. P., Steinhardt P. J., 2000, ApJ, 
530, 17 

Wetterich C, 1988, Nucl. Phys. B., 302, 668 

Wiaux Y., Vielva P., Martmez-Gonzalez E., Vandergheynst P., 

2006, Phys. Rev. Lett., 96, 151303 
Wiaux Y, Jacques L., Vandergheynst P., 2005, ApJ, 632, 15 
York D. et al. (SDSS collaboration), 2000, Astron. J., 120, 1579 



APPENDIX A: THEORETICAL DIRECTIONAL 
SPHERICAL WAVELET COVARIANCE 

When using azimuthally symmetric wavelets the spherical har- 
monic coefficients of the spherical wavelet transform may simply 
be written as the product of the signal and wavelet harmonic coef- 
ficients. This is not the case when using directional wavelets (i.e. 
wavelets that are not azimuthally symmetric). For this reason the 
theoretical form of the azimuthally symmetric wavelet covariance 
estimator derived by Vielva et al. (2006) is not trivially extended 
to directional wavelets. We derive here the theoretical directional 
spherical wavelet covariance. For simplicity, and without loss of 
generality, we ignore beam and pixel window functions in the anal- 
ysis (which, in any case, may easily be put in by hand at the end if 
required). 

Firstly, we derive an expression for the real space angular cor- 
relation function: 

oo f 

< n(co) f{w') > = Yj C ™ E Ye "' (a>) F ^ (fc/) 

;=o m=-e 

= f^-P e (co-co')Cf, (Al) 

where we have made use of (7) and also the addition theorem for 
spherical harmonics: 

J] Y, m (cj) Y* m (aj') = p £(aj . W ') 9 ( A2) 

m=-t 

where Y im (cu>) are the spherical harmonic functions and Pt{x) are 
the Legendre polynomials. 

Using (Al) we may write the theoretical wavelet covariance 

as 

< W"(a,b,p)W^(a,b,p)> 

I I dn(w)dO(w')C^Mfep( w ') 
x < n(a>) t*(o) f ) > 

= V^licff f dflMdfl(w') 

X KjbJV) taJbA *') P ^ ' W ') ' ( A3 ) 

Assuming ergodicity, we may relate (A3) directly to our wavelet 
covariance estimator. Furthermore, under our assumption that both 
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the CMB and galaxy density fields are isotropic, the wavelet covari- 
ance is independent of our choice of the position and orientation of 
the wavelet, hence we may set the rotation of the wavelet to zero, 
i.e. p = 0. Thus we may write 



Xf{a,b,y). 



00 21 + 1 C 

Y——C™ dfl(w)GV)«» (A4) 
hi 4n Js 2 



where 



G\oS) : 



X. 



(A5) 



Notice that the theoretical wavelet covariance specified by (A4) is 
independent of the orientation of the analysis y. The assumption of 
isotropy does not imply that the fields considered cannot contain lo- 
calised anisotropic structure, thus individual wavelet covariance es- 
timators obtained from the data will vary over orientations, but the 
theoretical prediction for each orientation will be the same when we 
consider the statistics of the aggregate fields. Essentially, we get a 
number of samples of a statistic for which we have a single the- 
oretical prediction, thereby actually increasing the performance of 
any subsequent comparison between the data and theoretical pre- 
dictions. Directional wavelets are therefore still advantageous in 
probing localised oriented structure in the data. 

Now we consider the spherical harmonic coefficients of G e (a>). 
Since Pt(oj • g/) is azimuthally symmetric it should be possible to 
write the spherical harmonic coefficients of G £ (oj) simply in terms 
of the harmonic coefficients of the wavelet: 



G\, m , = r 

Js 2 



(A6) 



(A10) 



where the last line follows from representing if/* ab (d>) by its har- 
monic expansion and noting again the orthogonality of the spher- 
ical harmonics (A8). For the case of azimuthally symmetric 
wavelets the m-modes of the spherical harmonic coefficients of 
the wavelet are non-zero only for m = 0, thus (A10) reduces to 
the form given by Vielva et al. (2006) (the (2£ + 1)/4tt discrepancy 
arises since Vielva et al. (2006) use Legendre, rather than spherical 
harmonic coefficients). We have therefore obtained a simple form, 
that may be computed easily, for the theoretical wavelet covariance 
of directional wavelets in terms of the wavelet spherical harmonic 
coefficients and the cross-power spectrum. 



where 



\ \ dQ(aj) 6Sl(o)') P e (cD • a>') Y* ffn/ (cj) F w (cd') 
Js 2 Js 2 



f f 

Js 2 Js 2 



dQ(aj) d£Ko>') 
Js 2 Js 2 21 + 1 

e 

X Y (m (w)Y* m (co')Y; m ,(co)Y t „ m „(co') 

m=—£ 

4tt 

Su'&U"5m'm" , 



2£+\ 



(Al) 



where we make use of the addition theorem for spherical harmonics 
again. The last line follows from the orthogonality of the spherical 
harmonics: 



dQ(cj) Y im (aj) Y € > m > (cj) = 6a>6 



Thus, the harmonic coefficients of G £ (co) are given by 



G £>m> ~ 



4tt 
2FTl 



(&ajb)e f m' $W 



(A8) 



(A9) 



We are now in a position to derive the final expression for 
the theoretical wavelet covariance in terms of the wavelet spheri- 
cal harmonic coefficients and the theoretical CMB -galaxy density 
cross-power spectrum. Expanding G € ((jj) in (A4) into its spherical 
harmonic decomposition yields 



/_n JS 2 



' €m (aj)i//* ab (aj) 



